Critical Behavior of Period Doubling in Coupled Area-Preserving 

Maps 

Sang-Yoon Kimtl 
Department of Physics 
•r^ ', Kangwon National University 

as _ 

On I Chunchon, Kangwon-Do 200-701, Korea 

> 

o 

:z; 



ON 






X 



Abstract 

We study the critical behavior of period doublings in N symmetrically cou- 



> 

(N 

O . 

(^ , pled area-preserving maps for many-coupled cases with A^ > 3. It is found that 



the critical scaling behaviors depend on the range of coupling interaction. In 
the extreme long-range case of global coupling, in which each area-preserving 
map is coupled to all the other area-preserving maps with equal strength, 



O I there exist an infinite number of bifurcation routes in the parameter plane. 



each of which ends at a critical point. The critical behaviors, which vary 
depending on the type of bifurcation routes, are the same as those for the 



H I previously-studied small A^ cases {N = 2,3), independently of N. However, 

for any other non-global coupling cases of shorter range couplings, there re- 
mains only one bifurcation route ending at the zero-coupling critical point, 
at which the N area-preserving maps become uncoupled, The critical behav- 
ior at the zero-coupling point is also the same as that for the small N cases 

{N = 2, 3), independently of the coupling range. 
PACS numbers: 05.45. +b, 03.20.+i, 05.70.Jk 
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I. INTRODUCTION 

Period doubling has been extensively studied in area-preserving maps of McMillan form 

xit + l) = -yit) + fixit)), 
y{t+l)=x{t), 

where {x{t), y{t)) is a state vector at a discrete time t and / is a nonlinear function. A well- 
known example is the Henon's quadratic map with f{x) = 1 — ax'^ [§]. As the nonlinearity 
parameter a increases, an initially stable orbit may lose its stability and give birth to a stable 
period-doubled orbit. An infinite sequence of such period-doubling bifurcations accumulates 
at a finite parameter value and exhibits a universal asymptotic behavior. However the 
asymptotic scaling behaviors for the area-preserving case [H-P] are different from those for 



the one-dimensional dissipative case |10 



An interesting question is whether the scaling results of area-preserving maps extend 
to higher-dimensional volume-preserving maps. Thus period doubling in four- dimensional 
volume-preserving maps has attracted much interest in recent years |^,|TT]-|r3| . It has been 
found in Refs. |r5|-P^ that the critical scaling behaviors of period doublings for two sym- 
metrically coupled area-preserving maps are much richer than those for the uncoupled area- 
preserving case. An infinite number of critical points form a critical Cantor set in the space 
of the nonlinearity and coupling parameters. The critical behaviors vary depending on the 
type of bifurcation routes to the critical points. These scaling results hold also for the 



three-coupled case [16 



It is also interesting to study whether or not the scaling behaviors for the two- and three- 
coupled cases extend to arbitrary A^-coupled (A^ = 2, 3, ■ ■ ■) cases. Thus we study the critical 
behaviors of period doublings in many coupled cases with A^ > 3 and compare them with 
those for the small N cases [N = 2, 3). In Sec. ^we introduce N symmetrically coupled area- 
preserving maps and discuss their general properties. Stability of orbits in this A-coupled 



map is also discussed in Sec. HI, We study the critical scaling behaviors of period doublings 



in Sec. |I^. It is found that the critical behaviors for many-coupled cases with A^ > 3 depend 
on the coupling range. In the extreme long-range case of global coupling, the stable regions 
of in-phase orbits of period 2" {n = 0, 1,2, ■ ■ ■) form a "bifurcation tree" in the space of 
the nonlinearity and coupling parameters, like the small A^ cases (A^ = 2,3). Consequently 
there exist an infinite number of bifurcation routes in the parameter plane, each of which 
ends at a critical point. The critical scaling behaviors, which vary depending on the type of 
bifurcation routes, are also the same as those for the small A^ cases (A^ = 2, 3), irrespectively 
of A^. However, for any other non-global coupling cases of shorter range couplings only one 
bifurcation route ending at the zero-coupling critical point is left in the parameter plane. 
The critical scaling at the zero-coupling critical point is also the same as that for the small 
A^ cases (A^ = 2,3), irrespectively of the coupling range. Finally, a summary is given in 
Sec. 

II. COUPLED AREA-PRESERVING MAPS 

In this section we discuss general properties of A^ coupled area-preserving maps in which 
the coupling extends to the Kth [1 < K < ^ {^^^) fo^ even (odd) A^] neighbor(s) with 
equal strength. 

Consider A^ symmetrically coupled area-preserving maps with a periodic boundary con- 
dition, 

^ [x(t + l) = -y(t) + F(x(t)), 
[y(t + l) = x(t), 

where x = (xi, . . . ,xi\f), y = (yi, . . . ,yN), F(x) = (Fi(x), . . . ,Fn(51)), and A^ is a positive 
integer larger than or equal to 2. Here Zm = {xm, Vm) ("^ = 1, ■ • ■ , ^) is the state vector 
of the mth element of the A^-coupled map T, and the periodic boundary condition imposes 
^m = Zm+N for ^H "^- The mth component Fm(x) of the vector- valued function F(x) is 
given by 



J \'^m) ~r QyXjm 3^m+l i • • • > ■^rra— 1 J; 



(2.2) 



where / is the nonhnear function of the uncoupled area-preserving map ( |1.1| ), and ^f is a 
couphng function obeying a condition 



g{x, . . . , x) = for any x. 
Thus the A^-coupled map T becomes: 

Xm{t + 1) = -ym{t) + f{Xm{t)) 

+ fi'(x„(t),...,x„_i(t)), 
ym(t + 1) = Xm(t), m = l,...,A^. 



(2.3) 



T: { 



(2.4) 



This A^-coupled map T is called a symmetric map |[13|-[16| since it has a cyclic permutation 
symmetry such that 



a ^Ta{z) = T(z) for all z, 



(2.5) 



where z = {zi, . . . , Zn), cr is a cyclic permutation of z such that crz = {z2, ■ ■ ■ ,zi), and 0"~^ 
is its inverse. The set of all fixed points of a is a two-dimensional (2D) subspace of the 2N 
dimensional state space, on which 



Xi 



xn, yi 



Vn- 



(2.6) 



It follows from Eq. (|2.5|) that the cyclic permutation a commutes with the symmetric map 
T, i.e., aT = Ta. Hence the 2D subspace becomes invariant under T, i.e., if a point z lies 
on the 2D subspace, then its image Tz also lies on it. An orbit is called an in-phase orbit if 
it lies on the 2D invariant subspace, i.e., it satisfies 



Xi{t) 

yiit) 



XNif) = X*{t), 

y^(t)=y*it) for allt. 



(2.7) 



Otherwise it is called an out-of-phase orbit. Here we study only in-phase orbits. They can 
be easily found from the uncoupled area-preserving map ( |1.1D since the coupling function g 
obeys the condition ( p.3|) . 



The Jacobian matrix DT of the A^-coupled map T is: 



DT 



DY 

In 



-In 




(2.8) 



where DF is the Jacobian matrix of the function F(x), In is the N x N identity matrix, and 
is the N X N null matrix. Since Det{DT) = 1, the map T is a 2N dimensional volume- 
preserving map. Furthermore, if DF is a symmetric matrix, i.e., DF* = DF {t denotes 
transpose), then the map T is a symplectic map because its Jacobian matrix satisfies the 



relation DT J DT = J 17 , where 



J 



V 







In 

N 



(2.9) 



The map T is also reversible since it can be factored into the product of two involutions 
1^, T = {TS)S [involutions satisfy (TS)^ = S^ = I (identity map)], where 



S: 



x(t + l) =x(t). 



(2.10) 



TS : 



(2.11) 



y(t + l) = -y(t) + F(x(t)), 

x(t + l)=y(t), 
y(t + l)=x(t). 

Such a decomposition represents a time-reversal symmetry because the map T is conjugate 
to its inverse T~^ by the involution S such that STS = T^^. The sets of the fixed points of 
the involutions, called the symmetry sets, are: 



Fix(^) = {(x,y)|y 



Ffxl 



}, Fix(T5) = {(x,y)|y = x}. 



(2.12) 



An orbit is called a symmetric orbit if it is invariant under the involution S; otherwise it 
is called an asymmetric orbit. It is easy to see that a symmetric periodic orbit has two 



symmetric points on the symmetry sets [T^: a symmetric orbit of even period has two 
points on one symmetry set and none on the other, but a symmetric orbit of odd period has 
one point on each symmetry set. Here we study only symmetric periodic orbits. 



Expressing the A^-coupled map ( p.4| ) in the form of second-order difference equations, we 
have: 



T : Xm{t + 1) + Xm{t - 1) 



Fm(x(t)) 
fiXmit)) 
+g{Xm(t),...,Xm-l{t)), 

m = l,...,N. 



(2.13) 



Consider an element, say the mth element, in this A^-coupled map T. Then the (m ± 6)th 
elements are called the Sth. neighbors of the mth element, where I < S < ^ ("^^y^) fo^ even 
(odd) A^. If coupling extends to the Kth neighbor (s), then the number K is called the range 
of the coupling interaction. 

A general form of coupling for odd A^ (A^ > 3) is given by 

K 



9[xi, 



.Xn 



2K + 



- ^ [u{xi+i) -u{xi)], 



l=-K 



K 



- XI u{xi+i) -u{xi) 



A^-1 



K = l,..., 



(2.14) 



where c is a coupling parameter and u is a, function of one variable. Here the coupling 
extends to the Ki\i neighbor(s) with equal coupling strength, and the function g satisfies 
the condition ( ^.3[ ) because the value of g for xi = ■ ■ ■ = xn becomes zero. The extreme 
long-range interaction for K = ^^^ is called a "global" coupling, for which the coupling 
function g becomes 



g{xi, • • • , a;Ar) = — X! Hxm) - u{xi)] 

m=l 

1 ^ 

— Y,u{x^) - u{xi) 

^^ m=l 



(2.15) 



This is a kind of mean-field coupling, in which each element is coupled to all the other 



N-l 



elements with equal coupling strength. All the other couplings with K < ^^^ (e.g., nearest- 

6 



neighbor coupling with K = 1) will be referred to as non-global couplings. The K = 1 case 
for A^ = 3 corresponds to both the global coupling and the nearest-neighbor coupling. 

We next consider the case of even N {N > 2). The form of coupling of Eq. (|2.14| ) holds 
for the cases of non-global coupling with K = 1, . . . , ^^-^ {N > 4). The global coupling for 
K = Y (^ ^ 2) also has the form of Eq. ( [^.15D , but it cannot have the form of Eq. ( [^.14[ ) 



since there exists only one farthest neighbor for K = ^^ unlike the odd A^ case. The K = 1 
case for N = 2 also corresponds to the nearest-neighbor coupling as well as to the global 
coupling, like the A^ = 3 case. 

III. STABILITY OF IN-PHASE ORBITS 

In this section we study stability of in-phase orbits in A^ coupled area-preserving maps. 
In case of global coupling the stability region of an in-phase orbit in the parameter plane 
is the same independently of A^, whereas for the other non-global coupling cases of shorter 
range couplings it depends on the coupling range K. 

The stability analysis of an orbit in many coupled maps is conveniently carried out 
by Fourier-transforming with respect to the discrete space {m} ||19|]. Consider an orbit 
{x(t)} = {xmii) ; m = 1, . . . , A^} of the A^ coupled map (|2.13|) . The discrete spatial Fourier 
transform of the orbit is: 

1 ^ 



A^ 



m=l 



j = 0,l,...,N-l. (3.1) 

The Fourier transform ^j(t) satisfies Cji^) = CN-j{t) (* denotes complex conjugate), and the 
wavelength of a mode with index j is ^ for J < y and j^ for j > y. 

To determine the stability of an in-phase orbit [xi{t) = ■■■ = XN{t) = x*{t) for all 
t], we consider an infinitesimal perturbation {6xm{t)} to the in-phase orbit, i.e., Xm{t) = 
x*{t) + 6xm{t) for m = 1,...,N. Linearizing the A^-coupled map ( ^.13 ) at the in-phase 
orbit, we obtain: 



where 



5Xm{t + 1) + Sx^{t - 1) = f'{x*{t)) SXmit) 



N 



+J2Gl{x*{t)) Sxi+m-lit), 



1=1 



(3.2) 



ax OXl+m-l 

dg{^) 



X\ = ---=Xiq=X 



dxi 

' Xl = --- = Xjsl=X 

Hereafter the functions G^'s will be called "reduced" coupling functions of g{x.) 
Let S^j{t) be the Fourier transform of 6xm{t), i.e., 

1 ^ 

^^ m=l 

j = 0,l,...,N-l. 



(3.3) 



(3.4) 



Then the Fourier transform of Eq. (|3.2|) becomes 



SC,{t + l) + 6Cj{t-l) 



N 

+5]G'Kx*(t))e2"('-^)^'/^] 
1=1 

j = 0,l,...,N-l. 



(3.5) 



This equation can be also put into the following form: 



/ 



50 (t + 
50 (t) 



1)^ 



/ 



L,(t) 



50 (t) 

^50(t-l) 
j = 0,l,...,iV-l, 






(3.6) 



where 



Hi) 



/ N \ 

' /'(x*(t)) + ^G'Kx*(t))e2-^('-iW/^ -1 * 



V 



(3.7) 



Note that the determinant of Lj is one, i.e., Det{Lj) = 1. 



Stabihty of an in-phase orbit of period q is determined by iterating Eq. (|3.6|) q times: 



y50(t + g-l) 



M, 



( \ 




^^ut--^)^ 


5 


J = 0,1,... 


,N-1 



(3.8) 



where 



t+q-l 

M, = n L,ik). 

k=t 



(3.9) 



That is, the stabihty of each mode with index j is determined by the 2x2 matrix Mj. 
Since Det{Mj) = 1, each matrix Mj has a reciprocal pair of eigenvalues, Xj and \J^. These 
eigenvalues are called the stability multipliers of the mode with index j. We also associate 
with a pair of multipliers {Xj, Xj^) a stability index, 

P^ = X. + Xt\ j = 0,1,...,N-1, 

which is just the trace of Mj, i.e., pj = Tr{Mj). 

It follows from the condition ( |2.3| ) that the reduced coupling functions satisfy 

N 



(3.10) 



J2Giix) = 0. 



(3.11) 



1=1 



Hence the matrix of Eq. (|3.71) for j = becomes: 



io(f) 



K 1 



A 







(3.12) 



/ 



This is just the Jacobian matrix of the uncoupled area-preserving map. Therefore the 
stability index po of the j = mode is just that for the case of the area-preserving map, i.e., 
Po depends only on the nonlinearity parameter a. While there is no coupling effect on po, 
coupling generally affects the other stability indices pj's (j ^ 0). 

In case of the global coupling of Eq. (|2.15| ), the reduced coupling functions become: 



GAx) 



;i - N)G{x) for / = 1, 
G{x) for / ^ 1, 



(3.13) 



where G{x) = -^u'{x). Substituting Gis into the (1, 1) entry of the matrix Lj{t), we have: 



N 

1=1 



2ni{l-l)j/N 



(3.14) 



for j = 0, 
—cu'{x) for j 7^ 0. 

Hence all stability indices p^-'s of modes with non-zero index j {j ^ 0) become real and the 
same, i.e., pi = ■ ■ ■ = Pn-i- Thus there exist only two independent real stability indices po 
and pi , the values of which are also independent of A^. 

We next consider the non-global coupling of the form ( p.l4| ) and define 



G{x) = 



2K + 1 



u{x), 



(3.15) 



where 1 < K < 



N-2 /A^-3^ 



G,(x) 



for even (odd) A^ larger than 3. Then we have 

-2KG{x) for I = 1, 

G{x) for 2 < / < 1 + K or 

for A^ + 1 - K < / < A^, 
otherwise. 



(3.16) 



Substituting the reduced coupling functions into the matrix Lj{t) of Eq. (|3.7]), the second 
term of the (1,1) entry of Lj{t) becomes: 



N 

^GKx)e2-(^-i)^/^ = -S^{K,j)cu'ix), 
1=1 



(3.17) 



where 



Sn{K,j] 



2K + 



4 ^ 

-L 7^ 1 



. 2^Jk . 

Sin -^-^ = 1 



k=l 



N 



g^n(2K + l)f 
{2K + l)sin^ 



(3.18) 



Hence all stability indices pj's with non-zero index {j ^ 0) become real, but they vary 
depending on the coupling range K as well as on the mode number j. Since Sn{K,j) = 
Siy{K, N — j), the real stability indices satisfy 



10 



Pj=PN-j, j = 0,l,...,iV-l. (3.19) 

Thus it is sufficient to consider only the case of < j < y i.^'^^) fo^^ even (odd) A^. 
Comparing the expression in Eq. ( p.l7[ ) with that in Eq. ( p.l4D for j 7^ 0, one can easily see 
that they are the same except for the factor Sn{K,j). Consequently, making a change of 
the coupling parameter c — > -g—n^, the stability index pj for the non-global coupling case 
of range K becomes the same as that for the global-coupling case. 

It follows from the reality of pj that the reciprocal pair of eigenvalues of Mj lies either on 
the unit circle, or on the real line in the complex plane, i.e., they are a complex conjugate pair 
on the unit circle, or a reciprocal pair of reals. Each mode with index j is stable if and only if 
the magnitude of its stability index pj is less than or equal to two {\pj\ < 2), i.e., its stability 
multipliers are a pair of complex conjugate numbers of modulus unity. A period-doubling 
(tangent) bifurcation occurs when the stability index pj decreases (increases) through —2 
(2), i.e., two eigenvalues coalesce at Aj = — 1 (1) and split along the negative (positive) real 
axis. 

When the stability index po for an in-phase orbit in an A^-coupled map decreases through 
—2, the in-phase orbit loses its stability via in-phase period-doubling bifurcation, giving rise 
to the birth of the period-doubled in-phase orbit. Here we are interested in such in-phase 
period-doubling bifurcations. Thus, for each mode with non-zero index j {j 7^ 0) we consider 
a region in the space of the nonlinearity and coupling parameters, in whcih both modes with 
indices and j are stable. This stable region is bounded by four bifurcation lines of both 
modes (i.e., those curves determined by the equations po = ±2 and pj = ±2), and it will be 
denoted by Un. 

In case of global coupling, those stable regions coincide, independently of A^ and j, 
because all stability indices pj's of modes with non-zero j are the same, irrespectively of A^; 
the stable region for this case will be denoted by Ug- An in-phase orbit is stable only when 
all its modes are stable. For the global-coupling case, Uq itself is the stability region of the 
in-phase orbit irrespectively of A^, because all modes are stable in the region Ug- 



11 



However the stable regions f/iv's vary depending on the couphng range K and the mode 
number j for the non-global coupling cases, i.e., Un = Un{K,j). To find the stability region 
of an in-phase orbit in an A^-coupled map with a given K, one may start with the stability 
region Uq for the global-coupling case. Rescahng the coupling parameter c by a scaling 
factor „ }^ ., for each non-zero j {j 7^ 0), the stable region Uq is transformed into a stable 
region Un{K,j). Then the stability region of the in-phase orbit is given by the intersection 
of all such stable regions [/at's. 

IV. CRITICAL BEHAVIOR OF PERIOD DOUBLINGS 

We are concerned about the critical scaling behavior of period doublings of in-phase 
orbits in A^ symmetrically coupled area-preserving maps ( p. 131) . Small A^ cases (A^ = 2,3) 
have been studied in Refs. |T^^T6|. It is interesting to study whether or not the critical 
scalings for the small A^ cases extend to the large A^ cases. Thus we study the critical 
behaviors in the large A^ cases and find that those for A^ > 3 depend on the coupling range 
K. For the A^-coupled map ( p^.l3D , we choose f{x) = 1 — ax"^ as the nonlinear function of the 
uncoupled area-preserving map and consider separately two kinds of couplings, the global 
and non-global coupling cases. 

A. Global Coupling 

We first study an A^-coupled map with global coupling, in which each element is coupled 



to all the other elements with equal strength. As shown in Sec. |IIT| , all stability indices 
Pj's of modes with non-zero j are not only the same (i.e., pi = ■■■ = Pn-i), but also 
independent of A^ for the global-coupling case. Thus the stability diagram of in-phase orbits 
in the space of the nonlinearity and coupling parameters becomes the same as that for the 
two-coupled case, independently of A^. That is, the stable regions of in-phase orbits of period 
2" (n = 0, 1,2, ■ ■ ■) form a "bifurcation tree" in the parameter plane. Consequently there 
exist an infinite number of bifurcation routes in the parameter plane, each of which ends at 

12 



a critical point. The critical behaviors, which vary depending on the type of the bifurcation 
routes, are also the same as those for the two-coupled case, irrespectively of A^ (for details 
of the iV = 2 case, refer to Refs. |T|-|l3). 

As an example, we consider a linearly-coupled case in which the coupling function ( |2.15| ) 
is 



g{xi,...,XN) = c 



1 ^ 

IT? Z^^m — Xi 



(4.1) 



Figure [| shows the stability diagram of in-phase orbits with period q = 1, 2, 4. These 
low-period orbits can be analytically obtained from the uncoupled second-order difference 
equation, x{t + 1) + x(t — 1) = f{x{t)) 0. 
The period-1 fixed point is 

1 



x*(0) = -(-1 + VT+a), (4.2) 

a 

and the independent stability indices are 

Po(a) = 2(l-v^rT^), (4.3) 

Pi{a,c) = poia)-c. (4.4) 

Its stable region in the parameter plane is limited by four bifurcation lines associated with 
tangent and period-doubling bifurcations of both modes (i.e., those curves determined by 
the equations pj = ±2 for j = 0, 1). 

When the stability index po of the fixed point decreases through —2, it loses its stability 
via in-phase period-doubling bifurcations, and gives rise to the birth of an in-phase orbit of 
period 2, 



x*{0) = -(1- Va-3), x*(l) = -(1 + Va-3). (4.5) 

a a 

Two independent stability indices of the period-2 orbit are 

Poia) = -4a + 14, (4.6) 

Pi(a,c) = po(a) + c(c + 4). (4.7) 

13 



Note that two period-2 stable branches bifurcate out of the period- 1 stable region in the 
parameter plane (see Fig. 0). 
An in-phase orbit of period 4, 

x*(0)=f^^)\ x*(l) = ^, 
V a\ a j \ a 



a;*(2) = -a;*(0), x*(3) = a;*(l), (4.8) 

bifurcates from the in-phase period-2 orbit at a = 4. Its two independent stability indices 
are 

Po(a) = 16a(2v^-a) + 2, (4.9) 

pi(a, c) = po(a) + c^ + ^\fa(? + 4(2v/a - \)(? 

-8v^(2a - 4yS + l)c. (4.10) 

As shown in Fig. |l|, each of the two period-2 stable branches also bifurcates into two period-4 
branches. Thus the stability region for this case consists of 4 stable branches. 

However it is difficult to obtain analytically higher periodic orbits with period larger 
than four. Thus one usually resorts to numerical computation. Successive bifurcations of 
stability regions have been also observed for the case of higher periods |T3|,|TB[. That is, each 
"mother" stable branch bifurcates into two "daughter" stable branches successively in the 
parameter plane; hereafter we call the direction of the left (right) branch of the two daughter 
branches L (R) direction. Consequently the stability region of the in-phase orbit of period 
2" (n = 0, 1, 2, . . .) consists of 2" branches. Each branch can be represented by its address 
[oq , . . . , a„] , which is a sequence of symbols L and R such that Qq = L and Oj = L or i? for 
i > 1. 

An infinite sequence of connected stable branches (with increasing period) is called a bi- 



furcation "route" |T3|,|T5|. Each bifurcation route is also represented by an infinite sequence 
of symbols L and R. Hence there exist an infinite number of bifurcation routes in the param- 
eter plane. A bifurcation "path" in a bifurcation route is formed by following a sequence 
of parameters {an,Cn), at which the in-phase orbit of level n (period 2") has some given 
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stability indices {po,pi) |T3|JT5[| . All bifurcation paths within a bifurcation route converge 
to an accumulation point {a*,c*). Here the value of a* is always the same as that of the 
accumulation point for the area-preserving case (i.e., a* = 4.136166 803 904...), whereas 
the value of c* varies depending on the bifurcation routes. Thus each bifurcation route ends 
at a critical point {a*,c*) in the parameter plane. 

The nonlinearity-parameter values a„ geometrically converge to a limit value a* in the 
limit of large n, i.e., 

a„ — a* ~ (5f " for large n, (4.11) 

where the value of Si is always the same as that of the scaling factor 6 (= 8.721 . . .) for the 

uncoupled area-preserving case, independently of bifurcation routes. However the scaling 

behaviors associated with the coupling parameter depend on the type of bifurcation routes. 

Consider p {p = 1,2, . . .) subsequences {cpm+j] m = 0, 1, 2, . . .} {j = 0, . . . ,p — 1) oi the 



coupling-parameter sequence {c„}. Then a bifurcation route is called a "period-p" route |TH 



if the p subsequences geometrically converge to a limit value c* in the limit of large m, i.e., 

Cpm+j — c* ~ 5^"^ for large m, (4-12) 

where the value of S2 is independent of j, but it depends on the type of bifurcation routes. 
As an example, consider the case of the lowest period-1 bifurcation routes [|1^-|15|. There 
exist three kinds of period-1 routes, called the S, A, and E routes. An S route is followed if 
one goes asymptotically in the same direction {L or R direction). Hence its address has the 
form, [u, (L, )°°] (= [u, L, L, . . .]) or [v, {R, )°°], where u and v are arbitrary finite sequences 
of L and R. The coupling-parameter scaling factor for this case is 62 = 4.000. Unlike the case 
of S routes, the direction of an A route asymptotically alternates between the L direction 
and the R direction. Hence its address has the form, [w, {L, R, )°°], where w is an arbitrary 
finite sequence of L and R. The value of 62 is —2.000 for all A routes except for the case 
of one special route, called the E route. The address of the E route is [{L, R, )°°] (i.e., the 
sequence w for this case is empty). (Hereafter only the routes with non-empty w sequences 
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will be called A routes.) The limit value of the coupling parameter for the E route is c* = 0, 
and the scaling factor is ^2 = —4.404. 

In order to see the periodicity of a bifurcation route, a "route sequence" has been in- 



troduced in Ref . |T^ . It is uniquely determined by the address of the route as follows. If 
the nth element of the address is the same as the next {n + l)th element, we assign a to 
the nth element of the route sequence; otherwise we assign a 1. Thus the route sequence 
becomes an infinite sequence of two symbols and 1. 

As an example we examine the asymptotic patterns of route sequences for the case 
of the period-1 routes. The route sequence of an S route is [u', (0,)°°] because one goes 
asymptotically in the same direction in an S route; here u' is an arbitrary finite sequence 
of and 1. On the other hand, the route sequence of an A or the E route is [ly', (1, )°°] 
because the direction of an A or the E route asymptotically alternates between the L and R 
directions. Here w' is an arbitrary non-empty sequence of and 1 for the case of A routes, 
whereas it is empty for the E route. Note that the route sequence of any period-1 route 
exhibits eventually a period-1 pattern. 

From the asymptotic period-p patterns of the route sequences for the low-period (p = 1, 2) 
routes, it has been conjectured in ||TB|] that if the route sequence of a bifurcation route exhibits 
an asymptotic period-p {p = 1,2,...) pattern, the route becomes a period-p one; otherwise 
it becomes a random one in which the coupling-parameter values c„ randomly converge to 
a limit value c* without any periodicity. 

B. Non-global Coupling 

Here we study the non-global coupling cases with K < ^ i^^^^) ^oi even (odd) A^ larger 
than three. Of the infinite kinds of bifurcation routes for the global coupling case, only the 
E route ending at the zero-coupling critical point (a*,0) is left as a bifurcation route in the 
parameter plane for all the non-global coupling cases. The critcal scaling behaviors near the 
zero-coupling critical point are also the same as those for the N = 2 and 3 cases of global 
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coupling. 

As an example we consider a linearly-coupled, nearest-neighbor coupling case with K = 1, 
in which the coupling function is 

f f (2^2 - xi) for A^ = 2, 
g{xi,...,XN) = < (4.13) 

\^{x2 + xn- 2x1) for iV > 3. 

Note that the two- and three-coupled maps for A = 2, 3 correspond to the global coupling 
case. The reduced coupling functions G/'s {I = 1, . . . ,N) of g defined in Eq. ( |3.3|) are given 
by 



G2{x) = ^ = G{x), Gi{x) = -G{x) for A = 2, 



G2{x) = Gn{x) = - = G{x), 

Gi{x) = ~2G{x), Gi{x)=0 {ly^l,2,N) for A > 3. 



Substituting G/'s into the (1,1) entry of the matrix Lj{t) of Eq. ( ^I7|) , we obtain: 



(4.14) 



(4.15) 



A o n ,^.,^^ -C^,! for A = 2, 3, 

£G,(x*(t))e2-('-i)^/^ = ' (4.16) 

'=1 [-Sn{1,j)c for N>3, 

where 5ji is the Kronecker delta (i.e., 6ji is for j = and 1 for j = 1), and 8^(1, j) = 
^sin^j^. Making a change of the coupling parameter c — > „ ,^^ .^ for each non-zero j, the 
stability index pj for A^ > 3 becomes the same as that for the N = 2 and 3 cases of global 
coupling. 

For each mode with non-zero index j, consider the stability region Un{1,j), in which both 
modes with indices and j are stable. The stability region Ug for the global-coupling case is 
independent of both A^ and j because all stability indices p/s of modes with non-zero j are 
the same, independently of A^. Rescaling the coupling parameter c with the scaling factor 
^^-7^, the stability region Ug for the N = 2 and 3 cases of global coupling is transformed 
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into the stablity region t/Ar(l, j) for A^ > 3. Since the factor Sn is dependent on the index 
j, the stabihty region Un{1,j) varies depending on the index j. Then the stabihty region of 
an in-phase orbit, in which all its modes are stable, is given by the intersection of all such 
stable regions Un's. 

As an example we consider the case A^ = 4. Figure ^ shows the stability regions of the 
2"-periodic (n = 1,2,3,4) orbits for this case. We first note that the scaling factor „ A ., 
has its minimum value | at j = 2. However [74(1, 2) itself cannot be the stability region 
of the in-phase orbit of level 1 (i.e., n = 1), because bifurcation curves of different modes 
with non-zero indices intersect one another. As shown in Fig. ^(a), the branch including a 
c = line segment remains unchanged, whereas the other branch becomes flattened by the 
bifurcation curve of the mode with j = 1. Due to the successive flattening with increasing 
level n [see Figs, ^(a) and (b)], of the infinite kinds of bifurcation routes for the global- 
coupling case remains only the E route ending at the zero-coupling point (a*, 0). Thus only 
the zero-coupling point is left as a critical point in the parameter plane. Note also that 
Figs. 0(a) and (b) nearly coincide near the zero-coupling point except for small numerical 
differences. 

We now examine the scaling behaviors near the zero-coupling critical point for the A^ = 4 
case of nearest-neighbor coupling. Consider a self-similar sequence of parameters {an,Cn), 
at which the in-phase orbit of level n has some given stability indices, in the E route for the 
global-coupling case. Rescaling the coupling parameter with the factor |, this sequence is 
transformed into a self-similar one for the A^ = 4 case of nearest-neighbor coupling. Thus 
the "width" of each stability region in the E route for the case of the global coupling is 
reduced to that for the A^ = 4 case of nearest-neighbor coupling by the scaling factor |, 



whereas the "heights" of all stable regions in the E route remain unchanged |2y]. It is 
therefore obvious that the scaling behaviors near the zero-coupling critical point for the 
nearest-neighbor coupling case are the same as those for the global coupling case. That is, 
the height and width /i„ and Wn of the stability region of level n geometrically contract in 
the limit of large n. 



/in ~ ^1 "■, Wn ~ ^2 " ^OT large n, (4.17) 

where 5i = 8.721 ... and ^2 = -4.404. 

The scahng results for the nearest-neighbor coupling case with K = 1 extend to all 
the other non-global coupling cases with 1 < K < ^ {^^) fo^ even (odd) N. For each 
non-global coupling case with K > 1, consider a mode with index jmax for which the factor 
Siy{K,j) of Eq. ( |3.18| ) becomes the largest one and its stability region Un{K,j^i^^) including 
a c = line segment. Here the value of jmax varies depending on the range K. Like the 
nearest-neighbor coupling case with K = 1, oi the two stable branches the one including 
the c = line segment remains unchanged, whereas the other one becomes flattened by the 
bifurcation curves of the other modes with non-zero indices. Thus the overall shape of the 
stability diagram in the parameter plane becomes essentially the same as that for the case 
with K = 1. Consequently only the i?-route ending at the zero-coupling point is left as a 
bifurcation route, and the scaling behaviors near the zero-coupling critical point are also the 
same as those for the global-coupling case. 

V. SUMMARY 

The critical behaviors of period doublings in N symmetrically coupled area-preserving 
maps are studied for many-coupled cases with A^ > 3. It is found that the critical scaling 
behaviors depend on the coupling range K. For the global-coupling case [K = ^ {^^) 
for even (odd) A^], the stable regions of in-phase orbits with period g = 2" (n = 0, 1, 2, . . .) 
form a "bifurcation tree" in the space of the nonlinearity and coupling parameters, like the 
small A^ cases (A^ = 2, 3). Hence there exist an infinite number of bifurcation routes in the 
parameter plane, each of which ends at a critical point. The critical behaviors, which vary 
depending on the type of the bifurcation routes, are also the same as those for the small A^ 
cases (A^ = 2, 3), independently of A^. However, for any other non-global coupling cases 
[1 < K < Y i^^Y^) for even (odd) A^], only the E'-route ending at the zero-coupling point is 
left as a bifurcation route in the parameter plane. The critical behaviors at the zero-coupling 
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critical point are also the same as those for the small A^ cases (A^ = 2, 3), independently of 
K. 
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FIGURES 
FIG. 1. Stability diagram of in-phase orbits in N linearly coupled maps with the global 

coupling. The horizontal (non-horizontal) solid and dashed lines denote the period-doubling and 

tangent bifurcation lines of the j = (1) mode, respectively. Note that all stability indices pj's of 

modes with non-zero j are the same (i.e., pi = ■ ■ ■ = Pn~i) for the global-coupling case. For other 

details see the text. 

FIG. 2. Stable regions of the in-phase period-2" orbits of level n = 1, 2, 3, 4 in four linearly 
coupled maps with the nearest neighbor coupling. The case of n = 1 and 2 (3 and 4) is shown in (a) 
[(b)]. Each stable region is bounded by its solid boundary curves. The period-doubling (tangent) 
bifurcation curve of the j mode of an in-phase orbit with period q is denoted by a symbol q^ . 
The values of the scaling factors used in (b) are 5i = 8.721 and 62 = —4.404. 
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